Mon. Not. R. Astron. Soc. OOO,ITH20lC2008l Printed 7 March 2013 (MN M£X style file v2.2) 



Relaxation and Stripping: The Evolution of Sizes, 
Dispersions and Dark Matter Fractions in Major and 
Minor Mergers of Elliptical Galaxies 

Michael Hilz 1,2 ,Thorsten Naab 1 *, Jeremiah P. Ostriker 3 ,Jens Thomas 4 , 
Andreas Burkert 2 and Roland Jesseit 5 

1 Max- Planck- Institut fur A strophysik, Karl-Schwarzs child- Str. 1, 85741 Garching, Germany 

2 Universitdts-Sternwarte Miinchen, Scheinerstr. 1 ,D-81679 Miinchen, Germany 

3 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544> USA 

4 Max- Planck- Institut fur extraterrestrische Physik, Giessenbachstrasse, D-85748 Garching, Germany 

5 Lifeware SA, Asylstrasse 112,CH-8032 Zurich, Switzerland 

Accepted ???. Received 111 in original form 111 

ABSTRACT 

We revisit collisionless major and minor mergers of spheroidal galaxies in the context 
of recent observational insights on compact massive early-type galaxies at high redshift 
and their rapid evolution on cosmological timescales. The simulations are performed 
as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure 
bulges as well as bulges embedded in dark matter halos. For major and minor mergers, 
respectively, wc identify and analyze two different processes, violent relaxation and 
stripping, leading to size evolution and a change of the dark matter fraction within 
the observable effective radius r c . Violent relaxation - which is the dominant mixing 
process for major mergers but less important for minor mergers - scatters relatively 
more dark matter particles than bulge particles to radii r < r c . Stripping in minor 
mergers assembles stellar satellite particles at large radii in halo dominated regions of 
the massive host. This strongly increases the size of the bulge into regions with higher 
dark matter fractions leaving the inner host structure almost unchanged. A factor of 
two mass increase by minor mergers increases the dark matter fraction within the 
effective radius by 80 per cent whereas relaxation in one equal-mass merger only leads 
to an increase of 25 percent. We present analytic corrections to simple one-component 
virial estimates for the evolution of the gravitational radii. These estimates are shown 
to underpredict the evolution of the effective radii for parabolic minor mergers of bulges 
embedded in massive dark matter halos. If such a two-component system grows by 
minor mergers alone its size growth, r c cx M a , reaches values of a w 2.4, significantly 
exceeding the simple theoretical limit of a = 2. For major mergers the sizes grow with 
a < 1. In addition, wc discuss the velocity dispersion evolution and velocity anisotropy 
profiles. Our results indicate that minor mergers of galaxies embedded in massive dark 
matter halos provide a potential mechanism for explaining the rapid size growth and 
the build-up of massive elliptical systems predicting significant dark matter fractions 
and radially biased velocity dispersions at large radii. 

Key words: galaxies: ellipticals - galaxies: evolution - galaxies: dynamics - galaxies: 
structure - methods: N-body simulations 



1 INTRODUCTION 

Recent observations have revealed a population of very 
compact, massive (« 1O 11 M0) and quiescent galaxies 
at z-~2 with sizes of about ~ lkpc (|Daddi et al.l 120051 : 
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Truiillo et all 120061: iLonshetti et al.l 


2007 


: iToft et al. 


2007 


; Zirm et al. 20071: Truiillo et al. 


2007 


; Zirm et al. 


2007 


; iBuitraEO et al.l 20081: Ivan Dokkum et alj 2008: 


Cimatti et all 120081; iFranx et all 12004 ISaracco et all 120091; 



Damjanov et al. 2009; Bczanson et al. 2009;| Newman et al.l 
20ld : ISzomoru et al.l [20121 : iMcLure et al.1 120121 '). Elliptical 
galaxies of a similar mass today are larger by a factor 
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of 4 - 5 (|van der Wei et all [2008) with at least an order 
of magnitude lower effective densities and significantly 
lower v elocity dispersions t han t he ir hig h- redshift counter- 



parts (Ivan der Wei et alj [iooEI. IgOoi 
20091 ; ICenarro fc Truullol 120091 : van Dokkimi et al 



2009 



Cappcllari et al 



van de Sande et al. l2011al : ISzomoru et al.l l201ll : 



Martinez-Manso et al.l 201 ll ). The measured small effective 



radii are most likely not caused by observational limitations, 
although the low density material i n the outer parts o f 
distant galaxies is difficult to detect (|Hopkins et al.l 12009). 
Their clustering, number densities and core properties 
indicate that they are probably the progenitors of the 
most massive ellipticals and Brightest Clu s ter G alaxies 
today l|Hopkins et al.l 120091 : iBezanson et al.l 120091 ). Such 
compact massive galaxies are extremely rare in t he nearby 



Valcntinuzzi ct al. 


2010; Raichoor et al. 


2012; Jiang et al.l 


2012; Fcrrc-Mateu et alj|2012l; iTruiillo et al.ll2012D. 



The high phase-space densities of these compact 
— 2 — 3 systems imply a formation mecha- 



dKeres et al. 20051; Naab et al.l 20071, 20091: Joung et al. 


2009 


Dekel et al. 


120091; Keres et al. 


2009 


1 Johansson et al. 


2009 


Oser et al.l 


20101; IWuvts et al 


2010 


; Bournaud et al. 


2011 


Oser et al.l 20121; Johansson et al.l |2012|). Frequent 



stellar minor and major mergers, as also expected in 
a cosmological context, are promising physical mecha- 
nism to explain the subsequent rapid size growth in 
the absence of significant add i tional dissipation and 



star formation (ICole et al.l I2OO0I; IKhochfar fc Silkl 120061; 



De Lucia et al.l 120061; iGuo fc White! 120081; IBezanson et al 



2009 1; [Hopkins et al 1120101 : ITruiillo et al.ll201ll ; ICooper et a! 



201 



2! 



McLure et al. 201 



12) 



Observations and theoretical studies provide evidence 
that massive early-type galaxies undergo on average about 
one major merger since redshift ~ 2 but significant! 



more minor mergers per unit time ([Bell ct al. 



Khochfar fc Silk 
Lotz et akl 



2011 



120061 ; iBell et al.1 120061 ; iGenel et al.l 



:antly 
2006; 



2008 



However, merger rates are notoriously 
difficult to measure and the uncertainties are still large, 
but assuming a canonical value of one major merger since 
z — 2, this would n ot be sufficient to explain the ob- 
served size evolution ([Bezanson et al. 1 120091 ; iTruiillo et al.1 
l201ll ; IWilliams et al.ll201ll) . In addition, major mergers are 
highly stochastic and some massive galaxies should have 
experienced no major merger at all, and would therefore 
still be compact today. Direct numerical simulations in a 
full cosmological context support the importance of nu- 
merous minor mergers for the assembly of massive galax- 
ies whose dissipative formation phase is followed by a sec- 
ond phase dominated by stellar accretion (predominantly 
minor mergers) onto the g a laxy, driving the size evolu- 
tion (IKhochfar fc Silkl 12009 ; iNaab et al.1 120091 ; lOser et al.1 
l2010l ; lLackner fc Gunnl2012l ). This theoretical finding is sup- 
ported by recent direct observational and circumstantial ev- 
idence that h as recently been present e d in support of mi- 
nor me rgers (Ivan Dokkum et al.l 120101: ITruiillo et al.l l20lJ : 
IWilliams et al.ll201ll ; ICooper et al.ll20121 )~ 

Minor mergers are particularly efficient in reducing the 
effective stellar densities, mildly reducing the velocity dis- 
persions, and rapidly incr easing the sizes , building up ex- 
tended stellar envelopes |Miller fc Smith! Il98d : IVillumsenl 



1983 




2007 




2010 





Farouki et akll 19831; INaab et al.ll2009l; iBournaud et al 



Bezanson et al.ll2009|; IHopkins et~ 2010l: 



_ _ ...... _ Oser et al 

20121 ; lLaporte et alj|2012l ; lLackner fc Ostrikerll201(]| ) 



However, there are some doubts as to whether minor mergers 
alone are suffi cient or whether additional physical pr o cesses 
are required dNaab et al.1 120061; INipoti et al.1 120091 . |2012| ; 
iNewman et al.ll2012l ; ICimatti et al. 20121 ). Major mergers of 
ellipticals will contri bute to mass grow th and will change 
kinematic properties (INaab et al.| [2006). However, their im- 
pact on the evolution of stellar den sities, velocity dis- 
persions and sizes is relatively weak JWhitel H978I . Il979l ; 

iBovlan-Kolchin et al.ll2005l ; INipoti et aljl2009h . 

Using t he virial theorem. ICole et al.l (|2000| ) , INaab et al.l 

( 2009) and IBezanson et al.1 {2009) presented a simple ana- 
lytical estimate of how sizes, densities and velocity disper- 
sions of one-component collisionless systems evolve during 
parabolic, purely stellar mergers with different mass ratios. 
According to this simplified model, the accretion of loosely 
bound material (minor mergers) results in a significantly 
stronger size inc rease than predicted for major mergers 
(INaab et al.ll2009T ). Using a similar approach IBezanson et all 
(2009) argued that eight successive mergers of mass ratio 
1:10 can lead to a size increase of ~ 5, which corresponds to 
the observed difference between old compact galaxies and 
massive ellipticals today. Of course, this is only valid for 
global system properties like the gravitational radii and total 
me an square speeds. T he simple virial e stima tes presented 
by INaab et all l|2009l ); IBezanson et al.1 (|2009l ) did not in- 
clude violent relaxation effects like mass loss, oc curring dur- 
ing th e encounter or non-homology. However, ICiotti et al.l 
(2007) presented a more general study including the effect 
of gas dissipation and non-homology. 

Early papers on the interactions of spheroidal galax- 
ies already discussed many of the aforementioned effects us- 
ing N-body s i mulat ions of one-component spherical systems. 
IWhitel l|l978l . ll979h . who made one of the first simulations 
of this kind, already found that relaxation effects are im- 
portant in equal-mass encounters and change the internal 
structure of the remnants. The centr al regions contract and 
diffuse envelopes build up (see also iMiller fc Smith! Il980l ; 
IVillumsenlll983l ; iFarouki et al.lll983l ). leading to a break in 
homology. Furthermore, equal-mass mergers lead to the re- 
distrib ution of particles eventua l ly red ucing population gra- 
dients (|Whitel Il980l : IVillumsenl Il983l ). This is different in 
unequal-mass mergers, which can enhance popu lation gra- 
dient s by depositing satellite stars at large radii (|Villumsenl 
1983). However, these early models did not investigate the 
influence of an extended and massive dark matter halo. 

Although dissipationless minor mergers theoretically 
lead to an increase in sizes and a decrease of velocity dis- 
persions , it is not clear whe ther this scenario works quanti- 
tatively. iM^otreFlil] (|2003l ). using two-component models - 
argued that dry major and minor mergers alone cannot be 
the main mechanism governing the late evolution of elliptical 
galaxies, because their si mulated merger remnan ts did not 
follow the Faber- Jackson (|Faber fc Jackson l ll976l ) and Kor- 
mendy relations (|Kormendvlll977r ). although they stayed on 
the fundamental plane. Still, dry major and minor mergers 
can bring compact early-type galaxies closer to the funda- 
mental plane but the size increase estimated from idealized 
models takin g into account the cosmological context migh t 
be too weak (|Nipoti et aill2009l . libTl ICimatti et al.ll2012h . 
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In this paper we revisit collisionless one- and two- 
component equal-mass and minor mergers focusing on the 
merger dynamics and the effect of a dark matter halo on 
the resulting evolution of global theoretical and observable 
galaxy properties. In section [2] we present the method used 
to construct equilibrium one- and two-component systems. 
The effects of violent relaxation and stripping are discussed 
in section [3] followed by a detailed discussion of analytic es- 
timates for the structural evolution based on the virial the- 
orem in section [4] In section [5] we present the evolution of 
observable galaxy properties. We summarize and discuss the 
results in section [6] 



2 NUMERICAL METHODS 
2.1 Galaxy models 

For the initial galaxy models we assum e spherical symm et- 
ric, isotropic Hernquist density profiles (|Hernquistlll990h for 
the luminous as well as the dark matter component, 



Pi(r) 



Mi 



2tt r(r + a^) 3 ' 



(1) 



where pi, Mi and a,i are the density, the mass and the scale 
length of the respective component i. The potential is 



Mr) 



GMj 
r + cii 



(2) 



with the gravitational constant G. 

On the one hand the projected Hernq uist profile is a 
reaso nable approximation of the R 1 ^ law l|de Vaucouleursl 
1 19481 ) for the luminous co mponent (its Sersic ind ex how- 



ever is closer to n ~ 2.6, see iNaab fc Truiiilcl 



200611. On the 



other hand it is a good representation of the Navarro et al.l 
l|l997f ) profile for the dark matter component. Therefore we 
consider the Hernquist density distribution a sufficiently re- 
alistic description for the luminous and dark matter distri- 
butions of a typical elliptical galaxy. 

For simplicity, we assume isotropy of the velocity dis- 
tribution to construct a stable initial configuration, and the 
bulge and halo component are in dynamical equilibrium. We 
compute the distribution function (D F) /, for each com- 
pone nt i, using Eddington's formula (Bi nnev fc Tremainel 

iH), 



ME) 



d 2 p l d$T 



V^TT 2 J$ =0 0(*T VE - $T ' 



(3) 



where pi is the density profile of component i, E is the rel- 
ative (positive) energy and $t is the total gravitational po- 
tential <J?t = $*(+$ dm)- Solving distribution functions, is 
in general more complicated than using Je ans equations, but 
resul ts in more stable initial conditions l|Kazantzidis et all 
|2004 ). The distribution function of a two component Hern- 
quist model ca n be compu ted analytically, even in the 
anisotropic case l|Ciottilll996h . Our model, however, also al- 
lows for more general density distributions and we calculate 
fi numerically. As in our approach we have no analytic ex- 
pression for p;(<&t) (see Eq. [3} we have to transform the 
integrand of Eq. [3] to be a function of radius r, 



d 2 p 



d® 2 . 



dr 



dr 2 




100.0 



Figure 1. Top panel: Radial density profile for a Hernquist dis- 
tribution (160k particles). The black solid line illustrates the 
initial profile (t=0) and the dashed-dotted line the final profile 
(t = T ~ 100 X tdym where t^yn ls the dynamical time at the 
spherical half- mass radius). The analytical profile is indicated by 
the red dashed line. Inside 10% of the scale length (vertical solid 
line) the system is affected (increase in dispersion and decrease in 
density) by two-body effects. They conduct heat into the initially 
cooler central regions. However, overall, the system is dynami- 
cally stable for at least 100 dynamical times. Bottom panel: Initial 
(solid line) and final (black dashed line) radial velocity dispersion 
profile. 



d«I>T 
dr 



d 2 $ T dpi 
dr 2 dr 



dr 



dr. 



(4) 



This procedure always results in an analytical ex pression 
for th e integrand, even for more general 7- profiles ()Dehnenl 
Il993h with different slopes for the density distribution. As 
a consequence, the integration limits also have to change, 
e.g. $(r) = gets r = 00 and $(r) = E has to be solved 
(numerically) for the radius r. 

Once we have computed the DF we can randomly sam- 
ple particles with radii smaller than a given cut-off radius 
and random velocities, which are smaller than the escape 
velocity. Then the particle configuration for the galaxy is 
established using the Neumann rejection method. 

The one component model is described by two parame- 
ters, the scale length a* and the total mass M». For the two 
component models, including dark matter, we additionally 
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introduce the dimensionless parameters [s, and j3 for the scale 
length of the halo adm = /3a* and its mass Mdm = (iM, . 



2.2 Model Parameters and Merger Orbits 

For the total dark matter to stellar mass ratio we assume 
/i = Mdm/M* = 10 and the ratio of the scale radii is 
P = adm /a* = 11, for all simulations with two-component 
models. We perform a set of simulations for two different 
scenarios; both adopt M* = a* = 1.0. In the major merger 
scenario we simulate mergers of initially identical, spheri- 
cally symmetric one- or two-component models on zero en- 
ergy orbits. The encounters have parabolic orbits with and 
without angular momentum (head-on). For higher merger 
generations we duplicate the merger remnant after reach- 
ing dynamical equilibrium at the center and merge them 
again on orbits with the same energy but different infall di- 
rections. In total, we simulate three generations of head-on 
equal-mass mergers, and two generations with angular mo- 
mentum (see also Table [TJ. 

In the minor merger scenario our simulation sequences 
start with an initial mass ratio of 1:10, i.e. Mh os t = 10M sa t 
and a stellar scale radius of the satellite of a*, aat = 1.0. 
This choice for the satellite's scale radius seems unrealis- 
tically large but there are no direct observations for low 
mass galaxies at high redshift and the sizes of less massive 
ellipticals at present converge towar ds an effective radius 
of r e ~ lkpc l|Misgeld fc HilkeJl201 j ). Therefore the satel- 
lite galaxies have the same size as the compact early-type 
hosts, although they are an order of magnitude less massive. 
For comparison, we simulated two sequences (with one- and 
two-component models), of head-on minor mergers, where 
the satellite's scale radii are half the host's scale radius 
a*, sat = 0.5, though the satellites li e on an extrapolation 
of the observed mass-size relation of I Williams et all (|2010l l 
at z = 2 (see Table [I}. The host galaxy for the next gener- 
ation is the virialized end product of the previous accretion 
event. This host is merged with a satellite identical to the 
first generation. The mass ratio for this merger is now 1:11. 
We repeat this procedure until the host galaxy has doubled 
its mass, i.e. 10 minor mergers. The final mass ratio of the 
merger is 1:19. Again we simulate one- and two-component 
mergers with zero (head-on) and non-zero angular momen- 
tum. As the mergers of the bulge+halo model with angular 
momentum are computationally expensive we only simulate 
6 generations. 

For all head-on mergers we separate the centers by 
a distance d and assign them a relative velocity u rc i = 
2\J GMh/d, where Mh is the total attracting mass of the 
host galaxy within the radius d. This velocity corresponds 
to an orbit with zero energy and zero angular momentum, 
i.e. the galaxies will have a zero relative velocity at infinite 
distance. The distance d is always large enough to obtain 
virialized remnants at the end of each generation. As the 
merger remnants after the first generation will not be spher- 
ical anymore, their mutual orientation is randomly assigned 
at the beginning of each new merger event. For the mergers 
with angular momentum we set the impact parameters to 
half of the spherical half- mass radius of the host's bulge and 
separate the galaxies for enough so that the initial overlap 
is very small. 



Run 


Gen. 


a*, sat 


M ub (%) 




Blho 


3 


1.0 


12.3 


12.3 


Blam 


2 


1.0 


15.0 


15.0 


HBlho 


3 


1.0 


10.1 


2.5 


HBlam 


2 


1.0 


8.8 


2.0 


BlOhod 


10 


1.0 


21.8 


21.8 


BlOamd 


10 


1.0 


20.9 


20.9 


BlOhoc 


10 


0.5 


21.7 


21.7 


HBlOhod 


10 


1.0 


35.6 


20.4 


HBlOamd 


6 


1.0 


19.5 


7.9 


HBlOhoc 


10 


0.5 


20.9 


7.2 



Table 1. This table gives the name of the hierarchy (1st col- 
umn), the number of generations (2nd), the initial scale radius of 
the satellite (3rd) , the amount of unbound mass of the total final 
remnant (4th) and the corresponding stellar mass loss (5th). The 
name can be explained as followed; B/HB: bulge or bulge+halo, 
1/10: major/minor merger, am/ho: orbit with/without angular 
momentum. In the case of the minor merger scenarios, c/d indi- 
cates wether we chose a compact or diffuse satellite. 

2.3 Simulations and Stability Tests 

All simulations were per formed with VINE IWctzstci n et al.l 
2009; iNelson et al.ll2009T ). an efficient, parallelized tree-code. 
We use a spline softening kernel with a softening length e = 
0.02 for all runs. In genera l, the softening length dep ends 
on the particle number (e.g. lMerrit3ll996MDehnenll2001J ) and 
we found e = 0.02 to result in stable models. For the major 
merger simulations the host galaxy consists of TV* = 1.6 x 10 5 
bulge particles for the one-component (bulge only) model 
and TV* = 2 x 10 4 for the two-component model, which has 
an additional halo represented by JVdm = 2 x 10 5 particles of 
the same mass. For the accretion scenario, the one- and two- 
component host galaxies both have TV* = 10 5 bulge particles 
and the latter has Nom = 10 6 halo particles. The satellites 
have ten times fewer particles for all components. 

In Fig. [1] we demonstrate the stability of the bulge only 
model with 160k particles by comparing the initial and fi- 
nal (100 dynamical times) density (top panel) and radial 
velocity dispersion (bottom panel) as a function of radius 
to the analytical values. In general the model is very stable 
for the full simulation time, except in the innermost parts, 
where two-body relaxation becomes important. At the high- 
est densities, inside 10% of the scale radius the relaxation 
time Wax of the model is very small (i rc iax ~ 50 • tdyn). 
Consequently, two-body encounters change the central par- 
ticle's energy and deplete the high density regions. Looking 
at the initial and final number of particles within 0.1a» , we 
find that half of the particles escape this region and go to 
lower binding energies. However, at larger radii the models 
are very stable with a very good agreement to the analytical 
solution. The radii enclosing 30, 50 and 80 per cent of the 
stellar mass stay perfectly constant. 

The stability of the two-component system is demon- 
strated in Fig. [2] Our initial model, constructed of two Hern- 
quist spheres, is again very stable over a long simulation pe- 
riod of more than 60 dynamical times (Figure 2). The den- 
sity and the velocity dispersion do not change significantly. 
Again the innermost regions are affected by two-body relax- 
ation. However, looking at the mass radii of the bulge and 
the halo we observe no big changes. The apparent contrac- 
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Figure 2. Top left panel: Radial density profile for the two-component realization with 1100k particles. The bulge to halo mass ratio 
is 1:10. The solid black lines illustrate the initial profiles (t=0) of the bulge and the halo and the dashed-dotted lines their final profiles 
(t > GOtdyn)- The analytic Hernquist profiles are indicated by the red dashed line. Bottom left panel: Radial velocity dispersion of the 
total system (blue solid: initially, blue dashed-dotted: final), the bulge and the halo separately. Inside 0.3a» the model is affected by 
two-body relaxation, but overall it is stable. Right panels: Time evolution of the radii enclosing 80%, 50%, 30% and 10% mass (black 
lines from top to bottom) of the bulge (top panel) and halo (bottom panel). The red lines show the effective radius of the bulge (upper 
panel) and the gravitational radius of the whole system (bottom panel). Except for the 10% radius of the bulge, which shows a slight 
increase, all mass radii stay constant over > 60 dynamical times. 



tion of the mass radii enclosing 80% or 50% is less than 5% 
for the halo and less than 2% for the bulge. The gravita- 
tional radius and the effective radius which we use in this 
paper stay constant. Therefore we conclude that the regions 
we are interested in are not affected by two-body relaxation 
and other numerical effects of the initial conditions, and the 
results of our merger simulations should be robust for our 
choice of force resolution and particle numbers. 



3 RELAXATION & STRIPPING 

Our set of simulations indicate that dissipationless mergers 
of spheroids are strongly affected by two dynamical pro- 
cesses, violent relaxation which dominates the evolution of 
major mergers and tidal stripping which is important for mi- 
nor mergers. In the following we will discuss both processes 
in more detail. 



3.1 Major Mergers 

For equal-mass mergers, the effect of violent relaxation is 
very strong and has a significant effect on the differential 
energy distribution of the remnants. The left panel of Fig. 
[3] indicates, that the initial narrow distribution (black line) 
becomes broader with each generation. Tightly bound parti- 
cles become even more bound and some weakly bound par- 
ticles gain enough energy to escape the galactic potential. 
The theoretical framework of violent re laxation is ve r y com - 
plex and, since the pioneering work of iLvnden-Belll (|l967h . 
there have been many approaches to obtain a viable theory, 
which can describe the final e quilibrium c onfiguration of a 
violently relaxing system (e.g. IShul Il978l : iNakamural l200ol : 
lArad fe Johansson! |2005i). In the following we first give a 
short introduction to the original approach of ILvnden-Belll 
(|l967t) . before we discuss our results with respect to anoth er 
slightly different approach of lSpergel fc Hernquistl |l992). 
During the approach and interaction of two collision- 
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Figure 3. Left panel: Energy distribution for the initial one-component host galaxy and three generations of head-on major mergers 
(Blho). The initial distribution (solid black line) is broadened during each merger towards lower and higher (cscapcrs) energies. Finally 
the most bound particles have about two times their initial binding energy. The mean energy of the total system stays constant (vertical 
dashed lines). Right panel: Energy distribution for all particles of the progenitor galaxy (black solid line, same as in the left panel) and the 
same particles after one equal-mass merger (blue solid line). Particles at high (red solid, —0.85 < e < —0.81), intermediate (green solid, 
—0.55 < e < —0.51) and low initial binding energies (light blue solid, —0.1 < e < —0.06) are highlighted and their energy distribution 
after the mergers are shown by the dotted distributions. By violent relaxation the two innermost bins are broadened by violent relaxation 
to Gaussian distributions, with a width of a = 0.14. Their mean energies are shifted to lower values (short- and dashed vertical lines). 
The most weakly bound particles arc shifted to higher energies resulting in a significant fraction of escapers. 



less systems the total gravitational potential <E> varies with 
time, whi ch leads to a, non-conservation of energy of si ngle 
particles l| L vnden- Belli 1 1 9671 : ISpereel fc HernauistJ[l992l ). 

- - -— (5) 

dt ~ dt' ( ' 

where e is the energy per unit mass. In accordance with the 
time dependent virial theorem, 



(6) 



where I is the moment of inertia tensor, the galaxy will 
convert its total potential energy V into kinetic energy T 
and back. In equilibrium, I — so T = —E, V = 2E, with 
E — T + V being the total energy. Away from equilibrium 
the total energy E is constant, but T and V will vibrate 
about these values, which will widen the differential energy 
distribution N(E), where N(E) gives the number N of stars 
within an energy interval E + dE. 

This evolution is illustrated in the left panel of Fig. 
[3] where the energy distribution N(e) broadens with each 
merger generation v e ry s imilar to the results reported in 
ISpereel fc HernquistJ l| 19921 ). Using the Ansatz, that vio- 
lent relaxation can be approximated by scattering effects 
of single particles, they found an analytic prediction for 
the final equilibrium configuration for the energy distribu- 
tion of an equal-mass mergers where the scattering prob- 
ability functions become Gaussian. In the right panel of 
Fig. [3] we select particles in three different bins with low 
(-0.06 < e < -0.1), intermediate (-0.55 < e < -0.51) 
and high (—0.85 < e < —0.81) initial binding energies. Af- 



ter the final merger, the initially narrow low energy bins 
are significantly broader and the two most bound bins can 
be well fitted with a Gaussian of width a — 0.14. Bins at 
lower energies does not become a Gaussian, eventually due 
to escaping particles. The mean energies of tightly bound 
particles (red vertical lines) are shifted to even higher bind- 
ing energies. The temporal evolution of the latter energy bin 
is detailed in the left panel of Fig. [4] Until t — 100 there is 
little evolution since the galaxy centers are not yet interact- 
ing. Then, at the first encounter, the mean binding energy 
increases strongly and the distribution broadens by a fac- 
tor of 3 (a = 0.02 — > 0.06). When two galaxies fly apart 
the mean energy nearly reaches its original value without 
further broadening (t = 130). During the second close en- 
counter, this scenario repeats, i.e. the energy distribution is 
shifted to higher binding energies accompanied by a strong 
broadening (t = 130 — > 140), before it oscillates back into a 
less bound state (t = 140 — > 150). Now the particles reside 
at slightly higher than initial binding energies. In the right 
panel of Fig.|3]we depict the evolution of the mean potentials 
of the three energy bins in Fig. [3] which oscillate strongly 
for the tightly bound particles (red, green lines) and the en- 
ergy shifts and broadenings are obviously correlated to the 
potential fluctuations. On the other hand these fluctuations 
vanish rapi dly, due to phase mixing resu lting in incomplete 
relaxation l|Lvnden-B cll 1967; Shu 1978). Violent relaxation 
offers new energy states by deepening the total gravitational 
potential, into which the most bound particles are scattered 
making them even more bound. At the same time a fraction 
of the weakly bound particles reach positive energies and 
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Figure 4. Left panel: Temporal evolution of the energy distribution of highly bound particles (red bin in Fig. [3}. The corresponding 
time evolution during the merger of the average potential of this energy bin (red line) and the two others presented in Fig. \3\ The most 
bound bin shows the strongest fluctuation during the first (2, t = 120) and second (4, t = 140) encounter. The initial distribution at 
t = 100 (0) is slightly broadened to cruo = 0.02 (1) by two-body relaxation. At the first close encounter t fti 120 (2) the potential rapidly 
gets deeper (see also right panel), the particles are shifted to higher binding energies and the energy distribution widens to <ri2o = 0.06. 
As the two galaxies fly away from each other, the potential increases and the particles go back to lower binding energies (3) without 
further broadening of the distribution. During the second close encounter, the process is repeated, i.e. the distribution broadens when 
the potential gets deeper (4) and moves back to higher energies. After t = 150 (5) the central regions show only negligible potential 
fluctuations and the particle distribution is slowly affected by two-body relaxation. During the merger, the innermost energy bin is 
broadened from a = 0.02 to a = 0.1 by violent relaxation. In isolation, this energy bin is only moderately affected by two-body relaxation 
over the same time interval (a = 0.02 — > 0.03). 



escape the potential of the remnant. Two-body relaxation 
effects for the initial system, evolved in isolation, are much 
weaker (e.g. a = 0.02 -> 0.03 between t = 110 -> 160). 

In the top panel of Fig.[5]we demonstrate, that the over- 
all energy distribution for major mergers of two-component 
models (solid lines) evolves similar to one-component mod- 
els, i.e. the tightly bound particles go to states with even 
higher binding energies and some weakly bound particles be- 
come unbound (positive energies). However, with each gen- 
eration, the number of dark matter particles in the highly 
bound regions increases more than the number of bulge par- 
ticles. This behavior can also be seen in the bottom panel of 
Fig. [5] The fraction of dark matter to bulge particles con- 
verges to unity in the central regions. As the energy of the 
50% most bound bulge particles, indicated by the small ver- 
tical lines at the top of this panel, hardly evolves the intrin- 
sic structure of the system has to change. For example there 
are only very few dark matter particles in the most bound 
regions initially, but a non- negligible amount occupies the 
most bound state of the remnant. This implies that the sys- 
tem undergoes a 'real' change in terms of its dark matter 
fraction. The fact, that more dark matter than bulge parti- 
cles wander to higher binding energies is illustrated in Fig. 
[6]for the first merger of HBlho. The amount of dark matter 
scattered to more bound states is larger than the luminous 
matter for energies e ^ —IT, which is the region where the 
halo starts to dominate the mass (see top panel Fig. [5]). 

This demonstrates that violent relaxation rearranges 



the distributions of dark and luminous matter in energy 
space, which yields a higher dark matter fraction at the cen- 
ter of the final system. The observable consequences will be 
discussed in more detail in section [5~2l 



3.2 Minor Mergers 

In Fig.[7|and the top panel of Fig.[8]we show the total energy 
distributions (solid lines) for a sequence of head-on, one- and 
two-component minor mergers, respectively. In both cases 
nearly all escaping particles originate from the satellites (red 
dashed-dotted line in both figures), which indicates that vi- 
olent relaxation only affects the in-falling material and has 
negligible effects on the distribution of the host galaxy par- 
ticles. As the satellites are less bound than the host galaxy 
we find a very high fraction of unbound mass for the final 
remnants. Furthermore, we can see that almost no accreted 
particles assemble in the central regions. The loss of bind- 
ing energy of the most bound particles in both scenarios 
is caused by two-body effects, which conduct heat into the 
central high density regions (see also Fig.[T]). Combining this 
shift with the effect, that most satellite particles assemble 
at low binding energies, results in an increase of the mean 
binding energies. For the bulge only mergers, the decrease 
of the mean energy can also be predicted analytically. 

The potential energy for a Hernquist sphere is 
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Figure 5. Top panel: Energy distribution for the two-component 
initial condition (black) and after the first (blue), second (green), 
and third (red) generation of equal-mass mergers. The solid lines 
depict the distribution for all particles, the dashed lines for the 
bulge particles and the dotted lines of the halo particles. The gen- 
eral evolution of the total system is similar to the one-component 
system (left panel, Fig. [3j) . The mean binding energy of the sys- 
tems (vertical dashed lines) stays constant. However, relatively 
more dark matter particles than bulge particles are scattered to 
low energies (or higher binding energies), thereby changing the 
shape of the dark matter energy distribution and increasing the 
central dark matter fraction (see also Fig. ll5l l. Bottom panel: The 
ratio of dark matter to bulge particles (or mass as we have equal- 
mass particles) increases with each merger generation for energies 
lower than —0.8 and decreases at higher energies (only bound 
particles are shown) . Finally, there are almost as many dark mat- 
ter than bulge particles in the innermost bins. The small vertical 
lines at the top of the panel indicate the energy including 50% 
of the most bound bulge particles, which stays constant after the 
first merger. 



l|HernauistJll99d ), 



W = 



GM 
6a 



(7) 



and according to the virial theorem the total energy of a 
system in equilibrium is 

Additionally we define the mass ratio of accreted to initial 
host mass rj = and the ratio of the according scale radii 
as (" = Assuming energy conservation for a zero energy 
orbit, the system's final energy is: 

„ . „ Ml Ml 

Ef = 




Figure 6. This panel indicates the energetic origin of bulge and 
halo particles after the first merger (black to blue line in Fig. [5} 
in terms of particle fractions per energy bin. The lowest energies 
e < —2.0 were not occupied in the initial system but all parti- 
cles are scattered down from higher energies during the merger 
(indicated by <— and intermediate blue for halo and intermediate 
red for bulge particles). Energies e > —2.0 were already initially 
populated with particles and some of them stay in their initial 
energy bin (indicated by J, and dark colors). The light blue/red 
regions show particles scattered up from lower to higher energies 
(indicated by — >). We predominantly observe downflow in regions 
where the energy distribution has a positive gradients, i.e. mostly 
at high binding energies (Fig. [5]l . This process is driven by violent 
relaxation and scatters more dark matter than bulge particles into 
the center increases (bottom panel of Fig.0. Upflow is seen for 
regions with negative gradients in the energy distribution, most 
prominently at the lowest energy bins leading to escapers. 
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Furthermore we calculate the mean final energy e, 
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(ii) 

(12) 
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12a; 



I2a a 



(9) 



where we used the fact that for equal mass particles the 
total number of particles is equivalent to the total mass M. 
In the case of an equal mass merger of two identical systems 
rj = £ = 1 and Ef — Ei, i.e. the mean energy of the system 
stays constant (see also Figs. [3] [5}. But for the numerical 
setup of the first minor merger (BlOhoc), where r/ = 0.1 and 
C = 1, the final mean energy is £/ = 0.927 • e<, in agreement 
with the simulations (Fig[7]|. 

Taking a closer look on the energy distribution of the 
bulge of the scenario HBlOhod (dashed lines, top panel Fig. 
[8} we can directly see, that most stellar particles accrete 
at energies e > —0.4, creating an overdensity of bulge parti- 
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Figure 7. Differential energy distribution for the initial one- 
component system (black) and two (purple), five (blue), eight 
(green), and 10 generations of 1:10 head-on mergers (BlOhoc). 
The red dashed-dotted line indicates the energy distribution of 
all accreted material, which shows that nearly all escapers come 
from the satellites and nearly no particles assemble at the center. 
In contrast to equal-mass mergers (left panel, Fig. |4j all bound 
particles become less bound as the mean binding energy of the 
systems (vertical dashed lines) decreases with each merger gener- 
ation. 



cles. Consequently the ratio of dark matter particles to bulge 
particles 'decreases' for e > —0.4 (bottom panel of Fig. [S]). 
While this ratio stays constant for all particles with e < —0.4 
the lowest binding energy of the 50% most bound particles 
moves to higher energies and the dark matter fraction in- 
creases (see also Fig |15|) . Altogether we can see, that the 
dark matter halos of the satellites are dissolved very rapidly 
in the deep potential well of the host galaxy and efficient 
tidal stripping leads to a build up of an stellar overdensity 
at low binding energies. 

Comparing the results of minor mergers with either 
compact or diffuse satellites, the overall evolution stays the 
same for one- and two-component models. But, as the com- 
pact satellites have a higher binding energy, tidal stripping is 
less efficient and the particles assemble closer to the center. 
However, the very innermost regions are still hardly affected. 



4 VIRIAL THEORY 



Going back to virial theory, Na ab et al.l ([2009) presented 
a simple prediction of how stellar systems evolve during a 
merger event, which will be extended later on and therefore 
is repeated briefly. Using the virial theorem and assuming 
energy conservation we can approximate the ratios of the 
initial to the final mean square speed (vf/ f ), gravitational 
ra dius r g i /f and density p,/f o f a merging system. According 
to iBinnev fc Tremaine l|2008l ) the total energy of a system 



10 5 r 



10 4 r 



10° 



10* 



10 1 r 





10° 




100.0 


JM 




~a 
n 

Z 


10.0 


s 

•0 


1.0 


z 


0.1 




-2.0 



■1.5 



-1.0 



-0.5 0.0 



0.5 



Figure 8. Top panel: Energy distribution for the initial two- 
component system HBlOhod (black) and two (purple), five (blue), 
eight (green), and 10 (red) generations of 1:10 mergers. The solid 
lines depict the distribution for all particles, the dashed lines for 
the bulge particles and the dotted lines of the halo particles. 
The red dashed-dotted line indicates the energy distribution of 
all satellite particles. For minor mergers violent relaxation only 
unbinds satellite particles and the energy distribution of the host 
is almost unaffected. In contrast to equal-mass mergers (Fig. 
bound particles become less bound, due to two-body effects. Sim- 
ilar to one-component minor mergers (Fig.[7Jl, the mean binding 
energy of the systems (vertical dashed lines) decreases with each 
merger generation. Bottom panel: In stark contrast to equal-mass 
mergers (Fig. the ratio of dark matter to bulge particles at en- 
ergies e < —0.4 is unaffected. At higher energies the fraction of 
(stripped) bulge particles increases. The short vertical lines in- 
dicate the energy including 50% of the most bound bulge par- 
ticles. In contrast to equal mass mergers this energy constantly 
increases, i.e. the bulge grows into the halo such that the enclosed 
dark mass increases. 



1 . 2 . 1 GMf 

= --MilvT) = -. 
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(14) 



where Ei and Mi are the system's initial total energy and 
mass. The gravitational radius is denned as 



GMf 



(15) 



Ki + Wi = 



■Ki = -Wi 



with the total potential energy Wi. Now we define E a , M a , 
r g , a and (t>„) as the energy, mass, gravitational radius and 
mean square speed of the accreted system. Additionally to 
the former defined mass ratio 77 = M a /Mi, e = 
is the dimensionless velocity fraction. By combining these 
assumptions with equation [14] we obtain 

Ml 
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Figure 9. Top panels: Evolution of the mean square speeds (left), the gravitational radii r g (middle) and the spherical density within 
r g (right) for the equal mass merger generations. The circles (triangles) represent the evolution of the one- (two-) component models 
with (red) and without (black) angular momentum. The black short dashed lines indicate the prediction of Eqs. 1161181 Considering 
the escaping particles Eqs. I2U23I yield the long dashed lines and the dashed-dotted lines for the one- and two- component scenarios, 
respectively. For better comparison, the light grey circles highlight the evolution of the minor merger scenario BlOamd. Bottom panels: 
The symbols show the evolution of the minor merger scenarios (see Table[TJ for the one- (circles) and two- (triangles) component models 
with (red) and without (black) a ngular mom e ntum . The blue circles depict the results, using compact satellite galaxies. The dashed lines 
are the idealized expectations of lNaab et al.l l|2009h and the dashed dotted lines are accounting for the escapers in the scenario BlOhod. 
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for the ratios of the final to initial mean square speed, grav- 
itational radius and density. In the very simple case of an 
equal mass merger of two identical systems rj = e = 1, r g 
gets doubled fEq. H7|l . {v 2 } stays constant (Eq. ll6[) and p de- 
creases by a factor of 4 (Eq. I18p . If we use equations 7-9 for 
a minor merger scenario, where {v 2 } « (v'f) and e << 1, 
the size of the final system can increase by a factor of ~ 4, 
as r oc M 2 . Additionally, the final velocity dispersion and 
density are reduced by a factor of 2 and 32, respectively. 
These changes are quantitatively in good agreement with 
observations. 

However, this simple analytic model suffers from a num- 
ber of limitations, apart from the restrictions to parabolic 
orbits and collisionless sy stems. The effect of violent relax- 
ation (|Lvnden-Bel]||l967t ) in the rapidly changing potential 



during the merger will scatter particles in energy space, 
making some more bound and unbinding others from the 
system, thus energy is not perfectly conserved (see Section 
[3]). Additionally, realistic spheroidal galaxies are composed 
of two collisionless components, dark and luminous matter 
with different spatial distributions, which are expected to 
react differently to a merger event (|Naab fc Truiillol 120061 : 
iHopkins et~aHl2009l ). 



4.1 Major Mergers 

Starting with the major mergers, the top left panel of Fig. 
[9] shows the evolution of the mean square speeds of the 
one- and two-component equal-mass mergers with the to- 
tal bound mass of the system. According to Eq. [16] the 
mean square speed should remain unchanged (dashed line), 
but obviously it increases with each generation. As a con- 
sequence, the growth of the total gravitational radius (top 
middle panel of Fig. [9]) and the density decrease (top right 



Relaxation and Stripping 11 



panel of Fig. [9j are weaker than expected . The same trend 
was reported bv lNipoti et al.l l|2003U2C)09l ) who correctly ar- 
gued that the simple analytical prediction is only valid for 
an idealized case without escaping particles. In the 4th col- 
umn of Table [T] we show, that the amount of unbound mass 
after the merger is not negligible and adds up to about 12, 
15, 10 and 9% per cent of the total mass for the scenarios 
Blho, Blam, HBlho and HBlam, respectively. The same 
effect can be seen for the two-component case, where two 
merger generations with angular momentum have nearly as 
much mass loss as three generations of the head-on coun- 
terparts (w 9% compared to « 10%). Looking at the last 
column of Table [1] we can see that nearly all escaping par- 
ticles for the two-component models are from the halo, as 
nearly no stellar mass gets lost (M* lU b < 3%). 

We can now revisit the analytic pre dictions and take th e 
effect of escapers into account (see also lNipoti et al1 l2003). 
The energy equation using the energy of the bound final 
system Ef and the energy of the escaping particles E esc is 



Ef + E esc = Ei + E a . 



(19) 



We assume that the escaping particles have essentially zero 
potential energy, so that 



E e 



+ ^M esc (v 2 sc ) 



(20) 



With a = M cac /Mi as the ratio of mass, lost in escapers to 



initial mass and /3 



c ) / (Vi) as the ratio of the mean 



square speed of the escapers and the initial system, we can 
now re- write equations 1 161 to [181 as 
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The long dashed (one-component models) and dashed- 
dotted (two-component models) lines in the top panels of 
Fig. [9] indicate that the updated analytic predictions are in 
good agreement with our simulation results. The deviations 
are less than a few per cent for both models. 

The situation becomes more complicated if we separate 
the velocities of the bulge and the halo component (Fig. 
UOp . The mean square speed of the bulge (green squares) 
increases more (finally > 50%) with respect to the total sys- 
tem (black squares), whereas the halo (blue squares) speed 
stays below the total. Here violent relaxation and dynami- 
cal friction lead to an energy transfer from the bulge to the 
halo (see section [3}, i.e. the fi nal bulge is more tightly bou nd 
than the initial one (see also iBovlan-Kolchin et al]|2005l for 
a discussion of the effect of different orbits). 

This effect can be estimated based on the ratio of dark 
and stellar matter. The total kinetic energy of the system is 



m*(v t ) + mj, 



i{Vdn 



) = mtot(v t ot)- 



(24) 



With m to t = m* + rn dm and introducing A(v 2 / drn ) = 
( v »/dm) ~ («tot) we obtain, 

m*A{v 2 ) + m dm A{v dm ) = (25) 
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Figure 10. The evolution of the mean square speeds of the 
total system (black squares), the halo (blue squares) and the bulge 
(green squares) with bound bulge mass for two-component models 
on head-on orbits (HBlho). The more rapid increase of the bulge 
velocities can be explained by an energy transfer from the halo 
to the bulge (green solid line, see Eq. [26}. At the center this 
effect is even more efficient, as the mean square speeds within the 
spherical half-mass radius (corresponding triangles) of the bulge 
increase more. The dotted green line indicates the expectation of 
Eg. 1261 for the central region. 



and the additional growth of the stellar mean square speeds 
is 



(26) 



If we now add A{v 2 } to the mean square speed of the galaxy 
(v 2 ot) we can consistently predict the bulge dispersion (green 
solid line in Fig. I10|l . 

As violent relaxation scatters some particles into states 
with higher binding energy, the central region contracts 
slightly relative to the total system growth (depicted by the 
gravitational radius). In Fig. [TT]we show the radii enclos- 
ing the 20, 50 and 80 % most bound particles normalized 
to the evolution of the gravitational radius. The inner re- 
gions expand less, the outer regions more than the gravita- 
tional radius. This effect is consistent with the ana l ysis i n 
section [3] and was already described in IWhitd l| 19781 . il979l ). 
It leads to a significant non -homology of the system (see also 
IBovlan-Kolchin et al.ll2005t ), as it evolves through successive 
mergers. 

On the other hand, this high central density leads to 
higher central velocities (triangles, Fig. II Of) . If we now add 
A{v 2 ) (Eq. I26p to (v 2 ot ) for the central regions we again get 
a very good prediction of the central bulge velocity (green 
dotted line Fig. HP) . 

4.2 Minor Mergers 

The bottom left panel of Fig. |9j indicates that the mean 
square speeds of all minor merger hierarchies decrease with 
increasing mass. In all scenarios with diffuse satellites (black 
and red symbols), the evolution is very close to the virial 
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Figure 11. Evolution of the radii enclosing the 20% (dashed), 
50% (solid) and 80% (dashed-dotted) most bound particles nor- 
malized to the evolution of the gravitational radius , r g (black: 
Blho; blue: Blam; green: HBlho; blue: HBlam). The inner re- 
gions expand less and the outer regions more than the gravita- 
tional radius. All ratios are normalized to the initial values. 



expectations of Eqs. 1161 T8l (dashed line), although the mass 
loss is significant, especially for the two-component models 
(red and black triangles). In Table [1] we can see that the 
fraction of escaping particles is up to 35% for HBlOhod and 
around 20% for all other minor merger scenarios. Further- 
more, regarding the 2C models, most of the escape fraction is 
due to the dark matter particles. Going back to the evolution 
of (v 2 ), we can see that the corrected prediction of Eq. 1211 
(dashed-dotted line), which includes the effect of mass loss, 
perfectly fits the results (e.g. scenario BlOhod). Using more 
compact satellites the final decrease of velocities (blue sym- 
bols) is much weaker, because they are more tightly bound. 
As they have half the scale radius of the diffuse satellites, 
their binding energies and velocities are two times higher 
which then doubles the velocity fraction e = {v 2 } / '(v 2 ) of 
Eqs. 1 161 181 and yields a smaller decrease. In combination 
with the occurring mass loss, this explains the different evo- 
lution of the mean square speeds. Nevertheless, in all sce- 
narios the final mean square speeds of the total systems are 
10 — 30% lower compared to their initial host galaxies. 

The evolution of the gravitational radii (bottom middle 
panel, Fig. [9} of the six hierarchies evolve according to the 
mean square speeds, which is not surprising as r g oc l/{v 2 ) 
(see Eq. ll4p . In detail, this means, that the hierarchies with a 
diffuse satellite show a size increase, which is consistent with 
the analytic predictions of Eq. [17] (dashed line). As the com- 
pact satellites are not able to efficiently decrease the mean 
square speeds, the gravitational radii grow only marginally. 
However, for all minor mergers the maximum size growth is 
around a factor ~ 2.4. For completeness, the bottom right 
panel illustrates, that the mean density within the gravita- 
tional radius evolves according to the gravitational radius 
(p oc rj 3 ), i.e. p decreases at maximum by 90% (model 
BlOamd, red circles). 



5 OBSERVABLES 

5.1 Size &: Dispersion Evolution 

Now we switch focus from theoretical galaxy properties, like 
the gravitational radii, to directly observable galaxy proper- 
ties. We compute the effective radius, r e , as the mean radius 
including half of the projected bound stellar mass along the 
three principal axes, the mean projected stellar velocity dis- 
persion, o~ e , within this radius and E e as the mean projected 
mass surface density within r e . 



5.1.1 Major Mergers 

The black triangles in the top left panel of Fig. [T2] indi- 
cate, that the effective velocity dispersion a 2 for the two- 
component mergers with head-on orbits (HBlho) evolves in 
a similar fashion to the central mean square speeds of the 
bulge (green triangles Fig. 1 1 Op . with a final value a factor 
~ 1.9 higher than initially. In the case of one-component 
models, the central region of the remnants just suffer from 
the contraction effect (Fig. 1 1 1 f) and a 2 is only slightly higher 
than the mean square speed of the total system (top left 
panel Fig. |5J). This effect is stronger for the scenario with 
angular momentum orbits. Here, o 2 increases more com- 
pared to the head-on case (top left panel Fig. I12|l and the 
escapers change the dispersions of both systems with an ad- 
ditional transfer of kinetic energy from the halo to the bulge 
particles, (see section [2). 

In contrast to the gravitational radii, the observable ef- 
fective radii of one- and two-component major mergers (top 
middle panel Fig. I12p follow the simple analytic predictions 
of Eqs. 1 161 181 i.e. the effective radius grows proportional to 
the mass, r e oc M~ 1 ' . At first this seems surprising as from 
the virial theorem r e oc M/a 2 and we would expect smaller 
radii as a 2 increases. Projection effects can be ruled out, as 
the spherical half-mass radii of the bound remnants evolve 
similar to r e and it cannot be an effect of d ark matte r alone. 

As discussed in iNipoti et all l|2003h and 
iBovlan-Kolchin et al.1 (|2005l ) we also find, that the 
systems change their internal structure with each merger 
generation. This non-homology effect can be quantified by 
the structure parameter c which connects the stellar mass 
of the systems to its observed size and velocity dispersion 

M,=c-Ig£. (27) 

We use Al* as the bound bulge mass of the remnants to get 
a comparable value for c of both merger hier archies (see also 
IPrugniel fc Simien|[l997l ; INipoti et al1l2009l ). A change in c 
indicates that the merger remnants do not have a self-similar 
structure. For the major mergers we find a continuous de- 
crease of c with each merger generation (top panel Fig. I13[) . 
The decrease is stronger (a factor of 1.8) for the models 
including dark matter (triangles) than for the bulge only 
models (factor of 1.5; circles), which can be explained by 
an increasing central dark matter fraction. In section |3~T1 we 
have already shown, t hat this is not just an effect of increas- 
ing radii (as argued bv lNipoti et al.ll2009h . but due to violent 
relaxation. Our resul ts just confirm the fin dings for major 
mergers presented in IBovlan-Kolchin et al.l ( 2005). The top 
right panel of Fig. [12] depicts the stellar surface densities, 
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Figure 12. Same as Fig. [9] but for 'observable' bulge properties like the effective line-of-sight velocity dispersions (left panels), the 
projected half-mass (effective) radii (middle panels) and the effective surface densities (right panels) versus the bound stellar mass 
normalized to the initial stellar mass. The top panels show the evolution of the major merger and the bottom panels of the minor merger 
hierarchies, respectively. The light grey triangles in the top panels are the results of the minor merger scenario HBlOhoc. The dashed 
lines indicate the simple analytic predictions (Fas, fT6l 1171 and I18|l neglecting cscapcrs and energy transfer. The dotted line in the top 
left panel accounts for both effects fits perfectly crj for scenario HBlho (see also Fig. |10t . 



which evolve according to r e and decrease by almost an or- 
der of magnitude for a stellar mass increase of a factor ~ 8. 
For the size increase with mass 

r e xM", (28) 

we find a = 0.8 — 1.0 for the first equal-mas s merger genera- 
tion, similar to lBovlan-Kolchin et alJ ((2005 ) for comparable 
merger orbits. The dispersion evolves as 

M» oc erf , (29) 

with j3 = 3.3 - 5.1. 

5.1.2 Minor Mergers 

In the bottom panels of Fig. [12] we show the evolution of <r e 
(left), r e (middle) and E e (right) for all minor mergers. The 
stellar dispersion a\ evolves only weakly and does not show 
the theoretically expected decrease, except for the two bulge- 
only scenarios with diffuse satellites (BlOamd, BlOhod). 

On the other hand the effective radii r e (middle panel) 
grow significantly and in particular for the two-component 



models the final size can be a factor of 4.5 larger for a mass 
increase of a about a factor of two. This is more than the 
maximum value from simple virial expectation (Eq. ll7p . For 
all minor merger scenarios of two-component models we get 
a mas s evolution of r e oc M~ 2 ' 4 , similar to iNipoti et al.l 
(|2012t ) for comparable orbits. This rapid increase is caused 
by stripped satellite particles which assemble at large radii 
(see section [3]). As an example we show the final stellar 
surface density distribution of model HBlOnod in Fig. 1141 
The host particles are concentrated at the center similar to 
their initial distribution. The accreted stellar particles have 
a much shallower surface density distribution and build up 
an envelope which dominates at radii r ^ r c . The struc- 
ture parameter evolves strongly and very similar for all two- 
component systems. Its evolution mainly reflects the rapid 
size growth (black, red and blue triangles in the bottom 
panel of Fig. I13|) . This strong evolution is also seen for the 
effective surface densities (bottom right panel of Fig. 1 14f) . 
which can decrease at maximum by an order of magnitude 
for a size increase of a factor of two. 

For bulge-only minor mergers scenarios the structural 
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Figure 14. Surface mass densities of bulge particles along the 
major axis, normalized to the final effective radius for head-on 
minor mergers (model HBlOhod) of two-component models. The 
grey solid line indicates the initial host stellar surface density, 
which almost remains unchanged (black dashed line). Most of 
the satellite stars (dotted line) assemble with a shallower surface 
density profile and dominate the outer parts of the total surface 
density of the remnant at r > r c (black solid line). The accreted 
stars are responsible for the rapid size growth. 



until the final coalescence, it suffers more from tidal strip- 
ping and builds up the most extended envelope of all bulge 
only models, which then results in the highest size growth 
(red circles, Fig. 1 1 2 j) . This implies, that the calculation of <r e 
also includes particles outside the innermost regions, where 
the velocities are lower and the velocity dispersion within 
the effective radius decreases. 



Figure 13. Evolution of the structure parameter c relating the 
stellar mass with the effective radius and velocity dispersion c = 
(r e a 2 )/(Mi,G) with bound stellar mass for 1:1 mergers (top panel) 
and 1:10 mergers (bottom panel). For all cases c decreases and 
homology is not preserved. In particular the structure of the two- 
component minor mergers (triangles) changes rapidly with added 
stellar mass. The colors and symbols in both panels are the same 
as in Fig. [9] 



changes depend on the initial conditions. Minor mergers 
with compact satellites (BlOhoc) evolve nearly self-similar 
with only mild changes in the structural parameter (bot- 
tom panel of Fig. 1 13 p and a weak evolution in a e , r e and 
E e (bottom panels of Fig. I12[) . The observable parameters 
evolve very similar to the theoretical ones (bottom panels of 
Fig. [9}. The other two bulge only minor merger models have 
weakly bound satellites, which already loose most of their 
material in the outer regions of the host galaxy and also 
build up an extended envelope with a corresponding change 
in structural properties (see also black and red circles in Fig. 
I13p . On the other hand, the development of an extended en- 
velope boosts the size growth of a system. As the sequence 
BlOamd with an angular momentum orbit needs more time 



5.2 Dark Matter Fractions 

As already discussed in section [3] mergers of bulge- halo sys- 
tems in general increase the central dark matter fractions. 
For major mergers this process is driven by violent relax- 
ation changing the composition of the remnants whereas in 
minor mergers structure is almost unchanged but the 'ruler' 
which is the size of the stellar system covers larger radii with 
initially higher dark matter contributions. The 'observable' 
quantity here is the dark matter fraction within the stel- 
lar half-mass radius which is presented in Fig. [15] (see also 
iBarnabe et~ai1l201ll ). For major mergers the dark matter 
fraction increases with every merger generation, probably 
saturating for later generations. A stellar mass increase of a 
factor two results in an increase of the dark matter fractions 
by ~ 25% for the first generation and ~ 18% and ~ 5% 
for the second and third generation, respectively. For minor 
mergers the increase in dark matter fraction is almost linear 
to the mass increase and significantly more dramatic than 
for major mergers. A mass increase of a factor two results 
in an increase of the dark matter fraction by ~ 80%, which 
is three times more than for major mergers, without any 
signature for saturation. This implies that merging, minor 
mergers in particular, is a significant driver for the increase 
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Figure 15. The dark matter fraction within the spherical half- 
mass radius of the bulge, f^ m = M^ m / (Af<jm + A^*)> versus the 
bound stellar mass of the two-component major (filled triangles) 
and minor (open triangles) merger remnants. The colors are the 
same as in Fig. [5] For major mergers, the dark matter fraction 
increases by 50% for a mass increase of a factor ~ 8 (25% for 
a mass increase of a factor ~ 2). Due to the rapid increase of 
the effective radii, the dark matter fraction of the minor mergers 
increases by almost 80% for a mass increase of less than a factor 
2. 



of dark matter in galaxies which had to be more baryon 
dominated in the past if their assembly history was domi- 
nated by stellar merging. We note that for minor mergers 
the increase in dark matter fraction seems to be roughly 
proportional to the mass growth. 



5.3 Velocity Dispersion Profiles 

One prediction of violent relaxation is that during a merger 
the system evolve to a state of higher entropy and the more 
complete the violent relaxation is the higher the entropy, 
eventually approaching the isothermal sphere. In reality, 
phase mixing damps the potential fluctuations of a merger 
rapidly, and violent relaxation is incomplete, which results 
in a final equilibrium distribution which does not reach the 
highest possible entropy state (see I White fc Naravanl (|l987l ) 
for a detailed discussion on maximum entropy states in 
galaxies) . 

In Fig. [I6]we show the radial velocity dispersion profiles 
for major mergers of one-component systems. The radial ve- 
locity dispersion of the fina l remnant ( red line) can be well 
fitted with a Jaffe profile l| Jaffe! 1 19851 ) over a large radial 
range. In the inner parts this profile rese mbles the singular 
isothermal sphere l|Tremaine et al.lll994l ). The correspond- 
ing changes of the structural parameter (see section [4} are 
also given in the same figure. In the case of two-component 
models (Fig. I17|l we find the same result, i.e. the total (solid 
line) a nd also the bulge profile (dot ted line) approach a Jaffe 
profile. ISpergel fc HernquistJ || 19921 ) find the same trend for 
only one generation of a head-on equal mass merger. This 
leads to the following tentative conclusion. The evolution of 
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Figure 16. Top: Final radial velocity dispersions of the one- 
component merger generations (Blho). The colors depict the dif- 
ferent generations. The dashed-dotted line indicates the velocity 
dispersion of a Jaffe profile, which has the same scale length and 
mass as the last remnant. The latter profile resembles the inner 
parts of the singular isothermal sphere which is a very good fit 
to our last merger remnant. As the structure parameter c de- 
creases with each generation homology is not preserved. Bottom: 
Anisotropy parameter /3 (eq. I30P against radius of three genera- 
tions of one-component equal-mass mergers. As f) > for higher 
generations the remnants become radially anisotropic over the 
whole radial range. The half-mass radius rso is the radius of the 
sphere, which includes half of the bound system mass. 



the merger remnants towards a more isothermal distribu- 
tion is not unexpected since the entropy cannot decrease in 
any collisionless process. However, by how much, or if at all, 
the entropy increases is hard to anticipate since, formally, 
there exists no maximum-entropy state for each merger gen- 
eration. The comparison with the Jaffe model shows that 
violent relaxation in the centre is very effective. From one 
merger generation to the next both the total mass and en- 
ergy change. This provides new phase-space states with ever 
higher binding energy that are populated by violent relax- 
ation. As a result, the nearly isothermal region of the rem- 
nant extends deeper into the centre after each generation. 

In the bottom panels of Figs. 1161 117l we can see that the 
initially isotropic remnants become radially anisotropic for 
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Figure 17. The top panel shows the same as in Fig. 1161 for the 
total remnants (solid lines) and the bound bulges (dotted lines) of 
the head-on two-component mergers (HBlho). The bottom panel 
depicts the anisotropy parameter for the bulge. Here, rso is the 
spherical half mass radii of the total (top) and stellar (bottom) 
bound remnants. If the bulge is embedded in a dark matter halo, 
/3 is slightly higher compared to the one-component case. 



most of the galaxy r > O.l rsp. Already after the first merger 
the anisotropy parameter |Binnev fc Tremainell2008l ) 



2o1 ' 



(30) 



becomes positive indicating radial anisotropy. This is not 
surprising, as we only use o rbits with very small or ze ro 
pericenter distances (see also iBovlan-Kolchin fc M a 2004). 

Next we focus on the evolution of the velocity disper- 
sion profiles (Figs. 1181 119|) of the minor mergers. From the 
previous analysis we have seen that violent relaxation has no 
significant influence on the central regions of the host galax- 
ies. This is also reflected in the weak evolution of the cen- 
tral radial velocity dispersion profiles oy(r) of the one- and 
two-component minor mergers (top panels of Figs. [181 119[) . 
In particular for the one-component scenario (e.g. BlOhoc in 
top panel of Fig.ll8p oy(r) keeps the initial Hernquist profile 
shape (black dashed line) over the whole radial range. Fur- 
thermore, as /3(r) ~ (bottom panel Fig. I18[) all remnants 



Figure 18. Top panel: The radial velocity dispersion for the 
head-on minor mergers of one-component models (BlOhoc) stays 
constant over most of the radial range. Only in the very cen- 
tral regions, it increases slightly with each generation. The black 
dashed line is the initial Hernquist profile and the red dashed- 
dotted line the velocity dispersion of all bound accreted particles. 
Bottom panel: For the whole bound remnant, the velocity distri- 
bution stays perfectly isotropic, as the anisotropy parameter /3 
stays zero. Looking at the accreted material (red dashed-dottcd 
line) , it gets radially anisotropic with increasing radius. In both 
panels the radius is normalized to the spherical half-mass radius 
of the bound system. 



(solid lines) stay isotropic. The accreted stars only show an 
increasing radial anisotropy f3(r) > with increasing radius. 
The aforementioned effects are similar for all one-component 
minor mergers, but slightly less pronounced in the cases with 
angular momentum. 

The merger remnants of two-component models (top 
panel Fig. I19|) also indicate characteristics which are con- 
sistent with the evolution of their differential energy distri- 
bution (Fig. [8]). The central regions are hardly affected and 
therefore, the total velocity dispersion profile (solid lines in 
the top panel of Fig. I19|) and the one of the halo stay con- 
stant. The accreted stellar particles create a bump in the 
energy distribution (dashed lines in Fig.[H|, which gets more 
and more prominent with each subsequent generation. These 
accreted particles induce an increasing stellar velocity dis- 
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Figure 19. Top panel: The radial velocity dispersion of the 
total system (solid lines) for the head-on minor mergers of two- 
component models (HBlOhod) stays constant over the whole ra- 
dial range. The dispersion of the bulge system (dotted line) builds 
up a prominent peak which comes from the accreted material, 
which is mostly affected by dynamical friction. The radii are nor- 
malized to the spherical half-mass radius of the bulge. Bottom 
panel: The anisotropy parameter of the bulge velocities gets ra- 
dially biased at radii greater than the spherical half-mass radii of 
the bulge. 



persion at radii larger than the spherical half-mass radius rso 
(dotted lines, top panel Fig. I19[l . As the final coalescence of 
the stellar component in the two-component scenario is on 
radial orbits, independent of the initial conditions (see also 
iGonzalez-Garcfa &: van Albadal ((2005)), the anisotropy pa- 
rameter becomes radially biased for all of our minor mergers 
(e.g. HBlOhod, bottom panel of Fig. [19]). This effect only oc- 
curs in the simulation including a dark matter halo, because 
then the angular momentum of the in-falling satellite is lost 
before the final merger, due to enhanced dynamical friction. 
Hence, most of the stellar particles approach on radial orbits 
and are, during the final coalescence, stripped before reach- 
ing the center. If we use the compact satellites, the overall 
trend does not change, but more material reaches the center, 
as the particles are more tightly bound and suffer less from 
tidal stripping. 



6 SUMMARY & DISCUSSION 

In this paper we have presented a series of numerical simu- 
lations of major (mass-ratio of 1:1) and minor (mass-ratio of 
1:10) mergers of spherical, isotropic early-type galaxies. The 
model galaxies consist of a bulge (one-component model) 
or a bulge embedded in a massive dark matter halo (two- 
component model) . After describing the procedure to set up 
the initial conditions we have demonstrated the stability of 
the one- and two-component galaxy models. By analysing 
the merger simulations we have identified two different pro- 
cesses dominating the evolution in sizes, velocity dispersion 
and dark matter fraction in major and minor mergers, re- 
spectively. 

Violent relaxation is the dominant mixing process in 
major mergers, broadening the energy distribution, unbind- 
ing weakly bound particles and making bound particles more 
bound. As a result the stellar velocity dispersion increases 
and the size increase is weaker (r tx M a with a < 1) than 
predicted fr om simple vi r ial ar guments, an effect already 
reported by iNipoti et al.l (|2003h . We present an analytical 
estimate for the size and dispersion evolution taking the 
mass-loss into account and also separating the effect on the 
dispersion of the stars and the dark matter halo for the 
two-component models. The physical origin of the increase 
of the dark matter fraction in equal-mass mergers reported 
bv lBovlan-Kolchin et al.l (120051 ) have been identified. In the 
transition region of luminous to dark matter domination, rel- 
atively more dark matter particles are scattered from higher 
to lower energies (becoming more tightly bound). This re- 
arrangement is driven by violent relaxation and leads to an 
increase of the the dark matter fraction within the observ- 
able effective radius of about 25 per cent for one equal-mass 
merger, equivalent to a stellar mass increase of a factor of 
two. 

Stripping of stars from the low mass satellites is the 
dominant process driving the size evolution in minor (1:10) 
mergers. The stripped stars are predominantly deposited at 
large radii, leaving the central structure of the host galaxies 
almost unchanged. This process was proposed to be the ori- 
gin o f abundance gradients in elliptical galaxies (|Villumsenl 
Il983t ) . The assembly of a 'halo' of accreted stars leads to a 
rapid increase of the effective radius. In systems with dark 
matter halos the stripped stars, with radially biased veloc- 
ity dispersion, are now populating regions where the host 
galaxy was previously dominated by dark matter. As a result 
the bulge size increase for two-component minor mergers is 
more rapid than predicted from simple analytical estimates 
(r oc M a with a > 2). This results in an about three times 
stronger increase of the dark matter fraction within the stel- 
lar half-mass radius (compared to equal-mass mergers) for 
a mass increases of a factor of two. 

Observational estimates of Ivan Dokkum et"aH (|2010h 
indicate a size growth of r e a M 2 04 , which would be qual- 
itatively consistent with two-component minor merger sce- 
narios. We have to note, however, that the exact details will 
depend on the structure of the galaxies, the merger orbits 
and the mass and extend of the dark matter halos etc., a pa- 



rameter spa ce too large to cover in this study. INipoti et al.l 
(2009, 2012) followed a similar approach in a cosmologically 
motivated context and found a much weaker size typical size 
increase (r e oc M 109 ). A possible reason for this discrepancy 
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is, that they averaged over all major and minor mergers. Ad- 
ditionally, they used a steeper central slope for the stellar 
density profiles leading to more bound satellites being even 
more compact than our satellites wh ich lie on an extrapola - 
tion of the z = 2 mass-size relation of lWilliams et al.l<|2010ft . 

As we, and others before, have shown, increasing dark 
matter fractions will change the ratio of dynamical to stel- 
lar mass. This effect, if relevant in nature, will contribute 
to the observed tilt of the fundamental plane, in addition 
to stellar populations and structural non-homology (e.g . 



Nipoti et al]|2003l; iBoylan-Kolchin et al.l 120051; iCiotti et al 
20071 ; iGraves fc Faberll2010l; iThomas et al.ll201ll ; iToft et a! 



2012ft . We also agree with iNipoti et al.l (120091 ) that the in- 
crease of the dark matter fraction is stronger for minor 
mergers and for this scenario is dominate d by t he rapid 
size growth. But in contrast to lNipoti et al.l |2009h . we find 
that the central dark matter fraction of equal-mass mergers 
represents a 'real' change in the internal structure caused 
by violent relaxation and not by the size growth of the 
galaxies. More massive galaxies, whic h are expected to have 



acquired more stars in late mergers toe Lucia et al 



2006 



2012 



Khochfar fc Silkl 120061 ; [Oser et aljfeoiol ; iMoster et al 
Laporte et alj|2012ft. will therefore hav e high er dark matter 
fractions (see also lLackner fc Ostrikerl (2010) for a detailed 
discussion). This picture is in agreement with recent high- 
resolu tion cosmological simulations (e.g. Ijohansson et al.l 
12012ft and might explain the potentially high dark mat- 
ter fractions measured for elliptical galaxies in particular 
if the merger history was dominated by minor mergers (e.g. 
Lackner fc Ostrikerl |2010| : iBarnabe etail l201ll ; I Auger et all 
2010ft . This explanation for the tilt of the fundamental 



plane is traditionally degenerate with uncertainties in the 
initial stellar mass functions. Some recent studies imply 
that the stellar mass-to-light ratio rises as a function of 
the stellar mass of the galaxies, leaving relatively little 
room for the presence of dark matter at the center of mas- 



sive galaxies (Ivan Dokkum fc Conroy||2010l; Cappellar 



Conrov & van Dokkum 20121; Ivan Dokkum fc Conro 



2012 



2012 



Ferreras et al.ll2012ft . If true this might be in considerable 



tension with cosmological simulations of massive galaxies 
and a merger driven size growth scenario - supported by 
massive dark matter halos - in general. 

The global mean square speeds of the systems de- 
crease in minor mergers but the 'observed' effective line-of- 
sight velocity dispersions hardly change. The most promis- 
ing realistic scenario investigated here, the two-component 
minor mergers, shows a relatively weak decrease in the 
effective velocity dispersion, eventuall y too weak to ex- 
plain the available observational data (ICenarro fc Truii Ho 
20091; Ivan Dokkum et ail 120091; Ivan de Sande et al.l l2011bF 



Toft et ail 2012T ) . However, reliable spectroscopic measure- 



ments are only available for few galaxies so far and the con- 
straints are weaker than for th e size evolution. In addition, 
studies like lCenarro fc Truiillol (|2009ft compared galaxies at 
a fixed stellar mass at different redshifts not taking the 
expected mass increase into account, i.e. massive galaxies 
at high redshift should be compared to even more massive 
galaxies (with higher dispersions) at present. In addition, 
details of the model, like mass concentration, orbits etc. will 
also affect the dispersion evolution. It is, however, reassur- 
ing that more realistic cosmological simulations, which also 



include a dissipative component, a re able to reprod uce the 
observed dispersion evolution (e.g. lOser et al.|[2012ft . 

In summary our work shows that dissipationless dry 
minor mergers of bulges embedded in dark matter halos 
can significantly increase the sizes of a compact spheroids 
resembling high-redshift systems. As the predictions for 
size growth are even above the observed values small 
amount of gas, which is known to reduce the size growth 
|Covington et al1l201ll ; [Hopkins et alj |2008). would perhaps 
not be enough to invalidate this scenario. Recent simulations 
of galaxy formation in a full cosmological context seem to 
support the conclusions drawn from the idealized simula- 
tions presented here with respect to the evolution of sizes, 
disper sions, abundance gradient as well as dark matter frac- 
tions ([Naab et alj|2007l. l2009[ ; Jjoung et al.ll2009l;IOser et al 
2010 1 |2012|; iFeldmann et all I2OI0I . l201ll ; Ijohansson et a! 
20121 ; lLackner et al.ll2012ft . 
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